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We reanalize the 77* — > 7r° transition form factor F- !T1 (Q 2 ) within the QCD light-cone sum rules 
method with the twist-4 accuracy. In computations, the pion leading twist distribution amplitude 
(DA) with two nonasymptotic terms and the renormalon method-inspired twist-4 DAs are used. 
The latters allow us to estimate impact of effects due to higher conformal spin components in the 
pion twist-4 DAs on the form factor Ft TJ (Q 2 ). Obtained theoretical predictions are employed to 
deduce constraints on the pion DAs from CELLO and CLEO data. 
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I. INTRODUCTION 



The form factor (FF) F 1T1 {Q 2 ) of the electromagnetic transition 77* — * tt° is one of the simplest exclusive processes, 
for investigation of which at large momentum transfers methods of the perturbative QCD (PQCD) can be applied 
[HQ- El- For computation of F ni (Q 2 ) various theoretical methods and schemes were proposed. They ra .2jp from the 
PQCD calculations 0,0], the light-cone sum rules (LCSR) |EIE0) to the running coupling approach ||, employed 
1— 1 ' to estimate power-suppressed corrections to F 7ri {Q 2 ). All of these methods are based on the PQCD factorization 
t-H , theorems, in accordance of which the amplitude of an exclusive process can be computed as the convolution integral 
of a hard-scattering amplitude and the process-independent distribution amplitude (DA) of an involved into the process 
hadron(s). The hard-scattering amplitude is calculable within the QCD perturbation theory, whereas the hadron (in 
■ our case, the pion) DAs are universal functions containing nonperturbative information on hadronic binding effects, 
<*S |j and cannot be obtained using tools of PQCD. Hence, the hadron DAs emerge as one of the two building blocks in 
JL 1 studying numerous exclusive processes. 

To calculate the 77* — * n° transition FF in the framework of the QCD LCSR method adopted also in this paper, 

1 the knowledge of the pion different twist DAs (p„(u,iip) is required. As we have just mentioned, they cannot be 
found by means of PQCD. Their factorization scale fjb F dependence is governed by PQCD, but an input information 
, at the starting point of evolution, i.e., the dependence of the DAs on the variable u (the longitudinal momentum 
fraction carried by the quark in the pion) at the normalization point /j,q, has to be extracted from experimental data 
or derived via nonperturbative methods, for example, the QCD sum rules [jj, instanton-based models |lCj . or lattice 
simulations Nevertheless, there exists the regular theoretical approach for treatment of the hadron DAs. It 

suggests the parametrization of the hadron DAs in terms of a partial wave expansion in conformal spin, and rely 
on the conformal symmetry of the QCD Lagrangian 12] . It is important that any parametrization of DA based 
on a truncated conformal expansion is consistent with the QCD equations of motion |T^ and is preserved by the 
QCD evolution to the leading logarithmic accuracy 0,0- Therefore, the conformal expansion provides a practical 
framework for modeling of the hadron DAs |!3lll4j and is widely used for investigation of numerous exclusive processes 
in QCD. 

Because of the increasing number of parameters at higher conformal spins and practical difficulties in phenomeno- 
logical applications, one has to restrict one's self by taking into account only the first few terms in the conformal 
expansion of DAs. As a result, the contributions of higher conformal spins to DAs in the existing calculations are 
neglected. At the same time, the suppression of higher spin contributions and the convergence of conformal expansion 
at present experimentally accessible energy regimes is not obvious and requires further studies. 

The renormalon model proposed in Refs. 0, 0] pursues to test precisely this issue; that is, to set a plausible 
upper bound for the possible contributions of higher conformal spins that so far escaped attention. The renormalon 
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approach employs the assumption that the infrared renormalon ambiguities in the leading twist coefficient functions 
should cancel the ultraviolet renormalon ambiguities in the matrix elements of twist-4 operators in a relevant operator 
product expansion. Such cancellation was proved by explicit calculations in the case of the simple exclusive amplitude 
involving pseudoscalar and vector mesons [Ti| . The idea of the renormalon model for the meson twist-4 DAs is to 
define them by taking the functional form of the corresponding ultraviolet renormalon ambiguities and replacing the 
overall normalization constant by a suitable nonperturbative parameter. It turns out that this is enough to obtain 
the set of two- and three-particle twist-4 DAs of the pion and and p-meson in terms of the corresponding leading 
twist (twist-2) DAs. It is remarkable that the set of twist-4 DAs, apart from the parameters in the leading twist 
DA, depend only on one new parameter. The latter can be related to the matrix element of some local operator and 
estimated using the QCD sum rules. 

A generic feature of the renormalon model is that it predicts higher twist distributions that are larger at the 
end points compared to the lowest conformal spin (i.e., the asymptotic) distributions, and are expected to modify 
a behavior of higher twist contributions in exclusive reactions. In fact, in our previous work [l7| we employed the 
renormalon-inspired twist-4 DAs for computation of the pion electromagnetic FF F 7T (Q 2 ) in the context of the LCSR 
method, and found that the new DAs enhance the twist-4 component of FF starting from Q 2 > 5.5 GeV 2 , and shift 
it towards larger values of Q 2 . Such modification affects the data fitting procedure and extraction of the parameters 
62(^0)' ^(Mo) m the pion leading twist DA, because in the renormalon approach the twist-4 contribution to F. K (Q 2 ) 
depends on the same parameters, as the twist-2 one, and is not a "frozen" background like in the standard analyses 



In the present work we reanalyze the 77* — > 7r° transition FF within the QCD LCSR method applying the 
renormalon-inspired model for the twist-4 DAs. We compare our predictions with the CELLO ^3 an d CLEO [2(j 
data on this process and deduce constraints on the parameters 62(^0)1 64(^0) a ^ tne normalization scale fi 2 , = 1 GeV 2 . 

This work is organized as follows: In Sec. II we define the two- and three-particle twist-4 DAs of the pion relevant 
to our present consideration and introduce their models in the renormalon approach. In Sec. Ill general expressions 
for the FF F 7T7 (Q 2 ) in the QCD LCSR method, as well as our results for the twist-4 contribution, are presented. 
In Sec. IV we compare our predictions with the CELLO and CLEO data on the 77* — * 7r° transition and obtain 
constraints on the parameters 62(^0), ^(Mo)' Section V contains our conclusions. Some important but cumbersome 
expressions are collected in the Appendix. 



In general, a pion is characterized by distributions of different partonic contents and twists. Its leading twist DA 
corresponds to a partonic configuration of the pion with a minimal number (quark-antiquark) of constituents. But the 
light-cone expansion of the relevant matrix element gives rise to two-particle higher twist DAs as well. The parton 
configurations with a nonminimal number of constituents (for example, quark-antiquark-gluon) are another source of 
the pion higher twist DAs. We concentrate here only on DAs that will be used later in our calculations. 

The light-cone two-particle DAs of the pion are defined through the light-cone expansion of the matrix element, 



Cl- 



ii. THE RENORMALON MODEL FOR THE PION DAS 



(tt°(p) 1^(2:2)7^75 [£2,2:1] u(xi)\0) 





(2.1) 



Here (f^(u, fi 2 F ) = ip^iu, (i F ) is the leading twist DA of the pion, and (p[ (u, p 2 ^), ^2(^1 Mf) are its two-particle 
twist-4 DAs. We use the notation [2:2, X\\ for the Wilson line connecting the points X\ and x 2 : 




(2.2) 



In Eqs. (|2.1|) and Ij2.2|) A = x\ — x 2 and u = 1 — u. 

The three-particle twist-4 DAs involving an extra gluon field can be introduced in the form |13| 

(t°(p) \u[-z) [-z,vz]j„~f 5 gG lip (vz) [vz, z] u(z)\ 0) 
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J Da . e -iP^-a 2 +a 3V ) ^(p^-p^)^^, 03,03) 



ZpPv 



ZpPv 



$x(ai, 0-2,0-3) f , 



pz / \ pz 

where the longitudinal momentum fraction of the gluon is 03 and the integration measure is defined as 

J Dai = J daida2<la^5{l — a% — a 2 — 03). 

The other pair of DAs is obtainable from Eq. after the replacement 7 5 G MP ->■ iG^ = le^^G^, 

7r°(p) u(-z) [-z,vz]'y l/ igGf J .p(vz) [vz,z]u(z) ; 



(2.3) 



(2.4) 
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(2.5) 



There exist one more three-particle twist-4 DA E^a-j) as well as, four quark twist-4 distributions, which we 
do not consider in this paper. 

The pion two- and three-particle twist-4 DAs are not independent functions, because the QCD equations of motion 
connect them with each other. From the analysis based on exact operator identities I2H]. it follows that 



(p^ (u) = I dv da\ 

Jo Jo Jo 



1-v 



da 2 — [2$_l - $11] (01,02,03), 
03 



ipf\u) + ^>{u) 
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uai — ua2 
1 a\ 



[2$ ± - $||] (oi,o 2 ,a 3 ), 



(2.6) 



where a% = 1 — ai — 02- 

The renormalon method provides the new, additional relations between different twist DAs of the pion. In order 
to explain principle points of the renormalon approach and derive relations between the pion twist two and four DAs 
in Ref. |Tfjj . the authors considered the gauge- invariant time-ordered product of two quark currents, 

(0 \T{d(x 2 )~(v"fh [X2,x\] 7r + (p)) , 

at small light-cone separations and expressed the matrix element in terms of two Lorcntz-invariant amplitudes 
Gi(u, A 2 ), £ = 1,2. They applied the operator product expansion to the amplitudes G;(it, A 2 ), computed the infrared 
renormalon ambiguities of the twist- 2 coefficient functions and ultraviolet renormalon ambiguities arising from higher 
twist operators, and proved that these ambiguities cancel exactly in OPE, rendering the structure functions Gi(u, A 2 ) 
unambiguous to the twist-4 accuracy. In the renormalon model, one defines the pion twist-4 DAs by keeping the func- 
tional form of the corresponding ultraviolet renormalon ambiguities and replacing the overall normalization constant 



plac 

cA 2 by the nonperturbative parameter S 2 /6. We refer the readers to Ref. [lg| for detailed analysis and calculations, 
and write down only final results: 



$±("1, oi 2 , o 3 ) = — 
6 
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$ll(ai,a2,a 3 ) = — 



(1-ai) 2 (l-a 2 ) 2 
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*j_(ai, a 2 , a 3 ) 



y7r(«i) ^(^2) 

1 — Qfj 1 — «2 



*ii(ai,a2,a3) = -^r 



a 2 <Ar(ai) . "1^(0:2) 



(1- Ql ) 2 (l-a 2 ) 2 



(2.7) 



In the left-hand sides of Eq. 1|2.7[) the substitution 013 = 1 — a± — a 2 is implied. 

Having substituted Eq. (|2.7() into Eq. (|2.6(l and computed the relevant integrals, one can obtain the pion two-particle 
twist-4 DAs tfii (u, i^p) and ip^\u, [i 2 F ) in terms of the leading twist DA. Stated differently, the twist-4 DAs (12.11) . 
Ij2.3(l . and (|2.5|l are determined solely by ip^{u, fi F ) and the new parameter 8 2 {^q). This parameter is related to the 
matrix element of the local operator 



5 2 {nl) ~ 0.2 GeV 2 



(2.8) 



and estimated from the 2-point QCD sum rules |22| . 

The last problem to be addressed here is a proper choice of the leading twist DA tp n (u, ftp). The renormalon 
method does not provide a prescription for that case, and we adopt a usual model for (pn(u, H-p) given by a truncated 
conformal expansion 

ip v (u, n%) = <p asy {u) [l+ J B 2 ( M |)C 2 3/2 (u-¥)+S 4 (M|)Cf 2 (u-u) + ---] , (2.9) 
and containing two nonasymptotic terms. In Eq. I|2.9[) <p a sy{u) is the pion asymptotic DA 

<Pasy (") = 6liU, 

and C?/ 2 (l;) are the Gegenbauer polynomials. The functions B n (fj, 2 7 ) determine the evolution of (f v (u, /i 2 ?) on the 
factorization scale /j, F , and at the next-to-leading order (NLO) are given by the following expressions: 



s 2 (/4) = + 



"s(Mf) , 2 /, 2 



4tt 



b A (v 2 )E 4 (v 2 F ) 



4tt 



where & 2 (/^§) and &4(/io) are the input parameters, which should be extracted from experimental data. 
In Eq. (TP E n ((j,%) and d k n {n 2 F ) (see Ref. 0) are defined as 



and 



"as(^F) 




"A) + /3iasQ4)/ 47r ~ 


(7i"//3i-7o7/3o)/2 


. "s(mo)_ 




./3o + /3ias(Mo)/47r_ 





7o * - 7o " - 2/3 



1 - 



qs040 

"s(Mo) 



(^-7o -2/3o)/2/3 ' 



(2.10) 



(2.11) 



(2.12) 



Here /3o, 7o and /3i, 7™ are the beta function and the anomalous dimensions one- and two-loop coefficients, respec- 
tively: 

2 38 
ft, = 11 --rc,, /3 1 = 102-yn / , 



n „ 2 100 4 728 

7o = 0, 7o = — , 7o = ^> 
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2 34450 830 4 662846 31132 

71 = °- 71 = "W " "sT 71 '' 71 = ~W " ^02T n/ ' (2 - 13) 

with rif being a number of active quark flavors. The functions d^if-tp) appear in NLO evolution formulas (|2.10|) due to 
mixing of partial waves in ip v (u, fi F ) corresponding to different conformal spins. The numerical values of the matrix 

M nk , 

M 02 = -11.2 + 1.73n/, M 04 = -1.41 + 0.565n/, M 24 = -22.0 + 1.65n/, 
and the standard two-loop expression for the QCD coupling, 

-in 



"s(a* ) = 



A)ln(M7A 2 ) 



In [ln(/i 2 /A 2 )] 
P 2 o ln( M 2 /A 2 ) 



(2.14) 



complete the necessary information on tp w (u, jJ> F ). 

The expansion of >p-n(u, up) in conformal spins (|2.9(l is the standard prescription in PQCD and is widely used in 
applications. However, to calculate the DAs tp^ (u, up), ip^ (u, Hp), as well as the twist-4 contribution to F^I^Q 2 ), 
we shall use the expansion of f v {u, jip) in powers of u, 

4 

¥>7r(u, = <Pasy{u) ^ K n{n F )u n . (2-15) 

n=0 

The coefficients K n (fi F ) in Eq. (|2.15|) are given by the following equalities: 

K a {p? F ) = 1 + 6B 2 (fip) + 15S 4 (a*f), KM) = -30 [B 2 (Mf) + 7B 4 (Mf)] , 

K 2 (pp) = 30 [B 2 (p 2 F ) + 28B4/4)] , #sG4) = -1260S 4 (/i|), 

# 4 (/40 = 63OB 4 (/40- (2-16) 

The distributions tp^ (u, pb F ) and i^ 4 "* ( u > Mf) m the renormalon approach were found in our work |17| using the 
expansion (|2.15|l 1 

4 

ip^ (u, n F ) = ^2K n (fi F )<pl(u), 

n=0 



^ 4 W20 = £tfn(/40*>nW- (2-17) 

n=0 

The explicit expressions of their components <fin(u) and </? 2 (it) are collected in the Appendix. 



III. THE FORM FACTOR F WJ (Q 2 ) IN THE QCD LCSR METHOD 

For the calculation of the electromagnetic 77* — * ir° form factor F 7T1 (Q 2 ) in the present work, we use the QCD 
LCSR method, which is one of the effective tools to estimate nonperturbative components of exclusive quantities 24]. 
The LCSR expression for the FF F^ 1 {Q 2 ) was derived in Ref. Q, where its tree- level twist- 2 and twist-4 components 
were found. The 0(as) correction to the twist- 2 part was computed in Ref. 0, and on the basis of these results 
constraints on the parameters b 2 (p^) and 6 4 (/^q) were extracted from CLEO data. This analysis was refined recently 
in Ref. where a new model for the pion leading twist DA was proposed. The renormalon approach to the twist-4 



1 The term — 61n 3 « 3 /45 i n ip \(u) (Ref. |17| . Eq. (2.18)) should be replaced by — 41m 3 m 3 /36. But this correction does not affect results 
and conclusions of Ref. Il 7l . 
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term leads to further insight on the FF, because it allows one to take into account effects due to the higher conformal 
spins neglected in the previous studies. 

The LCSR method is based on the analysis of the correlation function of the transition 7* (91)7* (92) — * t°(p) H 

dtxe-**' (n°(p) \T{j^x)M0)}\ 0) = ie^ a pq?q%F^(Q 2 ,q 2 ), (3.1) 

where Q 2 — —q 2 , q 2 = —q\ are the virtualities of the photons, j M = euHj^u + edd'y^d is the quark electromagnetic 
current and F nj » (Q 2 , q 2 ) is the form factor of the transition 7*(q , i)7*(q , 2) - > 71-0 (p)- 

For large values of Q 2 and q 2 , this correlator can be computed in PQCD. In the QCD sum rules method by matching 
between the dispersion relation in terms of contributions of hadronic states, which include a contribution of the low- 
lying physical states in the q-channel, i.e., that due to the vector p and to mesons, as well as a contribution coming from 
the continuum of hadronic states with the same quantum numbers, and the QCD calculation at Euclidean momenta, 
one can estimate the form factor i 7 ! n - 7 » (Q 2 ,q 2 )- 

After calculations, in the limit q 2 — > the formula for the FF Ft T1 {Q 2 ) can be obtained: 

F nj (Q 2 ) = — I ds ImF nj ,{Q*,s)ex V ( -f-^- I + - / — ImF ff7 . (Q 2 , s), (3.2) 

lrm p JO V M I 71 J so 3 




where 2 



-Im^ 7 .(Q 2 , S ) = 4^ 



(p.W(u,Q 2 ) 1 d<f>^(u,Q 2 ) 



Q 2 ds 



(3.3) 



In Eq. 1)3. 2(1 M 2 is the Borel parameter, m p is the mass of the p-meson (and p 2 = m% — 0), and So is the duality 
interval. By $( 4 )(it, Q 2 ) the following combination of the twist-4 DAs is denoted: 

^(u,Q 2 )=4(^(u,Q 2 )-^\u,Q 2 )) 



. , 1_u da 2 
aai I 

a 3 



1 — 2u + a± — 0L2 



a 3 



^u(ai,a 2 ,a 3 ) - *m (ai, a 2 , a 3 ) 



(3.4) 



CtR — 1 — Ot\— OL2 



The first term in the right-hand side of Eq. ((3.3(1 gives rise to the tree-level twist-2 component F$(Q 2 ) of the FF, 
whereas the second one generates the twist-4 contribution F${Q 2 ). In general, the F W7 (Q 2 ) has the following form: 

F^(Q 2 ) = F$(Q 2 ) + F^\Q 2 ) + F$(Q 2 ), (3.5) 

where F^ 2 1 ' as \Q 2 ) is the 0(as) correction to the twist-2 term. The formulas that determine F^ as \Q 2 ) are cumber- 
some and not written down here. Their explicit expressions can be found in Refs. 0,0). 

The twist-2 and -4 components of the FF can be rewritten in the form convenient for further analysis Q: 



(3.6) 



and 



Q 2 F${Q 2 ) 



V2U \ 1 f\ d^(u,Q 2 ) 



— / du - 



du 



■ exp 



+ 



uM 2 M 2 



Q 2 -' 



d<S>^(u,Q 2 ) 



du 



(3.7) 



where u = Q 2 /(s Q + Q 2 ). 



2 Starting from Eq. fM% = Q 2 . 
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0,0 0,2 0,4 0,6 0,8 1,0 

U 

FIG. 1: The components <!>„ (u) of $^(u, Q 2 ) 13.81 as functions of u. For comparison, the asymptotic DA I3.1U1 is also 
shown (the dashed line). 

Our main interest is the twist-4 term and its calculation using the renormalon-inspired twist-4 DAs. The DAs 
entering to Eq. I|3.4|) are known and after some calculations we get 



<P^(u,Q 2 ) = J2K n (Q 2 )<Pl?(u), 



(3.8) 



where the components $1 4 ' ) (u) are 



% ' (u) = 2S 2 [~2u 2 In u - 2u 2 In u] , $>\ 4) (u) = 2d 2 - [uu - 2u 3 In u - 2u 3 \nu] , 



^\u) = 25 2 



uu | 1 — q uu ) u2 (l — lnu + u 2 (l — 2u) In w 



$^ 4) ( W ) = 2S 2 \^ [36 - uu - I0u 2 u 2 ] + 2u 2 
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1- 2u- 



U 2 7i 3 
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In? 



2w 2 



1 - 2u + u 2 u 3 

5 



lni 



$i 4) (w) = 2<5 2 <^ uu 



4 7 _ 107 2 _ 2 

— I uu u u 

3 20 90 



20 9 

3 u + 5u 2 - 2u 3 

3 



In u + u 2 



20 9 o 

3 u + 5m 2 - 2m 3 

3 



In 7 



(3.9) 

The twist-4 function l|3.8|l coincides with one obtained in Ref. 25]. 

In Refs. the twist-4 contribution to F 7TJ (Q 2 ) was estimated using the asymptotic form of the relevant twist-4 

DAs, which lead to the simple expression 



2\ 2—2 



8 2 {Q 2 ) = 8 2 {ul) 



«s(Q 2 ) 



<*s(Mo) 



8C F /30 O 



(3.10) 



The function $ity(it, Q 2 ) and components of &( 4 \u, Q 2 ) are shown in Fig. ^ Without any detailed analysis, the 
difference in their behavior in the end-point regions u = 0; 1 is evident: we have emphasized that the renormalon model 
predicts DAs that are larger at the end points. The difference between them becomes more essential, when considering 
the behavior of d<&^ {u,Q 2 )/ du in the end-point regions. Thus, for the asymptotic model d^a) y {u, Q 2 )/du — at 
u = 0; 1. In the case of the renormalon-inspired DAs for n = — 4 we find d$n\u)/du = 4<5 2 (Q 2 ) at u — and 
d<S> ( n\u)/du = -4<5 2 (Q 2 ) at u = 1. 

The expression of the twist-4 contribution l|3.7|l through the spectral density p^(Q 2 , s) and that of p^(Q 2 , s) itself 
are presented in the Appendix. 
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0,1- 



b 2 °=0.2,b 4 °=-0.1 



0,0- 



0.8 



1,0 



1.2 



1.4 



M 2 (GeV 2 ) 



FIG. 2: The dependence of the scaled FF Q 2 F^(Q 2 ) on the Borel parameter M 2 . The DA with 62(^0) = 0.2, 64(^0) 
is used. For the solid curve Q 2 = 2GeV 2 , for the dashed curve Q 2 = 4GeV 2 , and for the dot-dashed one Q 2 = 8GeV 2 



-0.1 



IV. EXTRACTING THE PION DAS FROM THE EXPEIMENTAL DATA 



The LCSR expression for the pion electromagnetic transition FF can be used to extract constraints on the input 
parameters ^2(^0) ano - ^4(^0) °f the- leading twist DA. In order to perform numerical computations, we fix various 
parameters appearing in the relevant expressions. Namely, we take the Borel parameter M within the interval 
0.7 GeV 2 < M 2 < 1.4 GeV 2 and for the factorization and renormalization scales accept 

2 2 n2 

For the QCD coupling as(Q 2 ) the two-loop expression (|2.14|) with A3 = 0.34 GeV is used. The duality parameter 
so = 1.5 GeV 2 is determined from the two-point sum rules in the p-meson channel 26] . The normalization scale is 
set equal to ^ = 1 GeV 2 . We also use m p = 0.775 GeV and = 0.132 GeV. 

The Borel parameter dependence of the LCSR for different values of Q 2 is depicted in Fig. [3 In calculations the 
leading twist DA with &2(Mo) = ^2 = 0-2, ^4(^0) = ^4 = ~ 0.1, as well as the twist-4 function l|3.8|l are used. From 
this figure, one can conclude that the prediction for the FF is rather stable in the exploring range of M 2 . In what 
follows we choose the Borel parameter equal to M 2 = 1 GeV 2 . 

The scaled FF Q 2 F^ 7 (Q 2 ) and its different components are plotted in Fig. [3] As is seen, the leading order prediction 
for the FF in the QCD LCSR method considerably overshoots the data points. By including into consideration the 
NLO correction, one may only soften this discrepancy, but not remove it. In fact, the LCSR prediction with LO+NLO 
accuracy, for the values of the input parameters b®, 64 shown in the figure, again overestimates the experimental data. 
It is evident that DAs with 64 > will lead to a more great deviation from the experimental data than DAs with 
64 < 0. Therefore, the traditional treatment of the transition FF 0,0 would call for DAs with b\ < 0, because the 
asymptotic twist-4 contribution is not strong enough to compensate the growth of the LO+NLO term if 6° > 0. In 
the renormalon approach the situation, in general, remains the same. This fact is connected with the dependence of 
the twist-4 term on the input parameters 6°, 64. 

To understand this important point, it is instructive to explore the twist-4 contributions corresponding to different 
DAs. The relevant results are shown in Fig. 21 Here the dashed and solid curves are computed using the standard 
asymptotic function l|3.10|l and the renormalon-generated one (|3.8|) with 6" = ^4 = 0, respectively. The difference 
between them is evident: almost in the whole region of the explored momentum transfers Q 2 > 1.8 GeV 2 the 
higher conformal spin (renormalon) effects enhance the absolute value of the twist-4 contribution. The twist-4 terms 
corresponding to 6° 7^ 0, t|j 7^ (the dot-dashed and dot-dot-dashed curves) in the region 1 GeV 2 < Q 2 < 1.8 GeV 2 
are larger than the asymptotic contribution, whereas for Q 2 > 1.8 GeV 2 they run below both the solid and dashed 
curves. It is worth noting that at fixed &° = 0.2, curves b° 4 = ±0.1 (these values have been chosen as sample ones ) 
are close to each other: the sizeable difference between them appears only at Q 2 > 8 GeV 2 . This means that in the 
low-Q 2 region, DAs with b® = ±0.1 reduce the LO±NLO contribution almost in the same manner, and this effect is 
smaller than the corresponding effect in the case of the asymptotic contribution. In the domain of the high Q 2 , the 
twist-4 terms cut the LO±NLO result more effectively than the asymptotic twist-4 term, but because the LO+NLO 
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2 4 6 8 10 

Q 2 (GeV 2 ) 

FIG. 4: The twist-4 term as a function of Q 2 . The dashed curve is computed using the function &^} y (u, Q 2 ). The predictions 
for the twist-4 term obtained employing the renormalon-inspired DAs are shown by the solid and two broken lines. The 
correspondence between the lines and the input parameters is: the solid line, 62(^0) — 0, 64(^0) = 0; the dot-dashed line, 
b 2 {no) = 0.2, 64(^0) = -0.1; the dot-dot-dashed line, &2O0) = 0.2, 64 (/4q) = 0.1. 



(64 < 0) contribution itself is smaller than the LO+NLO (64 > 0) one, it terns out that only DAs with b\ < lead to 
agreement with the data. 

In general, calculations of Ref. 00] correspond essentially to the "minimal" model of the twist-4 effects, where 
the restriction to the lowest conformal spin probably underestimates the effect, while the renormalon model is a 
"maximal" model, where these effects are probably somewhat overestimated. Therefore, the renormalon model for 
the twist-4 DAs allows us to put a theoretically justified bound on the twist-4 contribution to the pion transition form 
factor. The change in absolute value of the twist-4 term is not dramatic, as it may be expected. To quantify this 
statement we introduce the ratio R(Q 2 ), 



R(Q 2 ) 



[F${Q 2 )] ien 



[F^ ) (Q2)]stand' 

and demonstrate its numerical results in Fig.[3]for some selected values of b®, b\. 
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8 10 

Q 2 (GeV 2 ) 

FIG. 5: The ratio R(Q 2 ) as a function of Q 2 . The solid line corresponds to the input parameters 62(^0) — ^(Mo) — 0- 
The dashed line describes the same ratio, but for 62(^0) = 0-2, ^(Mo) = 0.1, while the dot-dashed one corresponds to 
fe(Mo) = 0-2, 6 4 (^) = -0.1. 

In Fig. we plot the ltr area in the 6° — 64 plane, extracted at the scale fi^ = 1 GeV 2 in the result of the fitting 
procedure. One can see that this area stretches in the lower-half plane occupying a large region. The parameter 63 
varies within the limits 

b% £ [0.09, 0.37], 

whereas at some fixed b% the variation of b\ is even larger. For example, at b\ = 0.2, 64 takes values in the region 

b% e [-0.05,- 0.45]. 



0,0 


-0,2- 



-0,4- 




-0,6 

0,10 0,15 0,20 0,25 0,30 0,35 

K 

FIG. 6: The la area in the 6° — b% plane extracted from comparison of the CELLO and CLEO data with the LCSR prediction. 
The central solid rectangle denotes the point b% = 0.27, b% — —0.3. 
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FIG. 7: The form factor Q 2 F^(Q 2 ) as a function of Q 2 . The shaded area demonstrates la region for the transition FF. In the 
fitting procedure the solid data points have been used. The central solid curve corresponds to the parameters 6° = 0.27, b% — 
-0.3. 



The la region for the transition FF is shown in Fig. [7] In the fitting procedure we employ the CLEO data, and 
the two CELLO data points at Q 2 = 1.26 and 1.7 GeV 2 . The CELLO data have an important impact on the fitting 
procedure. 

The central values for the parameters bl an d bl extracted in this work are 

GeV 2 ) = 0.27, b° 4 (l GeV 2 ) = -0.3. (4.1) 



It is difficult to present the la region of the parameters by a simple formula, therefore we refrain from such attemps. 

As is seen the curve in Fig. 0in the low-Q 2 region decreases rapidly, deviating considerably from the CELLO point 
at Q 2 = 0.94 GeV 2 : the DA gdj describes the CELLO and CLEO data located only in the domain 1.26 GeV 2 < 
Q 2 < 7.9 GeV 2 . 

Evolved to the normalization scale /z 2 , = 5.76 GeV 2 , the parameters (|4.1|l become equal to 

b° 2 ~ 0.184, bl ~ -0.179. (4.2) 
Equation l|4.2|l can be compared with the predictions = 5.76 GeV 2 ) from Ref. |g, 

b a 2 = 0.19, b a A = -0.14, (4.3) 

and Ref. 

b% = 0.2, bl = -0.17. (4.4) 

From this analysis it becomes evident that both the standard and renormalon approaches predict the pion leading 
twist DA with a negative coefficient bl < 0. 

In our previous paper [17J we obtained the constraints on the pion leading twist DA from LCSR analysis of the 
pion electromagnetic FF. In calculations we employed the renormalon model for the twist-4 distributions (|2.17|) . but 
took into account the evolution of the pion DAs with the lower than in the present work accuracy (i.e., with the LO 
accuracy). At the scale ^ = 1 GeV 2 we found 

b% = 0.2 ± 0.03, bl = -0.03 ± 0.06. (4.5) 

The la area depicted in Fig. [6] overlaps with (|4.5() (it is not shown in the figure) in the region determined by the 
following values of the parameters 



bl ~ 0.2 - 0.23, b\ ~ -0.05 - (-0.09). 
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In general, an la area for two quantities is not necessarily equal to the overlap region of two independent la areas. 
Nevertheless, it is important, at least as the first approximation, to determine the pion DA that satisfactorily describes 
the both of these form factors. In order to get more precise estimate, F 7T (Q 2 ) has to be calculated with higher accuracy, 
and a joint treatment of F^ 1 (Q 2 ) and F^^Q 2 ) must be performed. These tasks are beyond the scope of the present 
work. 



V. CONCLUSIONS 



In this work we have used the pion renormalon-inspired higher twist DAs to calculate the twist-4 contribution to 
the transition form factor F^^Q 2 ). The renormalon method has allowed us to express two- and three-particle twist-4 
DAs in terms of the pion leading twist DA and one additional parameter 8 2 . In other words, in this method the 
twist-4 distributions are determined by the twist-2 one unambiguously, that restricts a freedom in the choice of DAs, 
increasing, at the same time, the predictive power and reliability of QCD results. 

The higher twist distributions introduced in Ref. [T^ | embrace higher conformal spin effects which so far escaped 
attention. They are larger at the end points and, as expected, give rise to larger higher twist effects in exclusive 
reactions. Our results on the twist-4 contribution to the transition FF F 7r7 (Q 2 ) confirm this assumption. Indeed, 
in the region of high momentum transfers the absolute value of the renormalon-generated twist-4 term exceeds the 
asymptotic one by a factor 1.5 — 2. Nevertheless, the LCSR results obtained using the renormalon-inspired and 
standard DAs predict the pion leading twist DAs with b\ < 0. Stated differently, effects due to higher conformal spin 
components of the twist-4 DAs remain under control and do not spoil the QCD LCSR method. 

The new contribution of this work is that the renormalon approach has allowed one to put an upper bound on the 
twist-4 contribution to the light-cone sum rules and obtain estimates of the effects due to higher conformal spins. On 
the example of the 77* — » ir° transition, it also demonstrated limits of theoretical uncertainties inherent to our 
knowledge of exclusive processes. 
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APPENDIX 

The components ip L (u ) and <Pn(u) of the pion two-particle twist-4 DAs ip^ (u, ip^ (u, /J, 2 ?) are given by the 
following expressions |l7j : 

ifio(u) — 5 2 |u [Inu — Li2(u)] + u [Inu — Li2(u)] — uu+ j ' 



(Pi(u) = S 2 <u 







u 




{■ 


[(- 


' 2 


4) 



Inu — Li2(u) 



u 
2 



Inu — Li2(u) 



o ± 2 2 

—uu -\ — u u H 

6 2 6 



ip\{u) = S 2 I u 



1 + u — —u ) Inu — Li2(u) 



2 

u — — u ) Inu — Li2(u) 



£ o 2 2 

-- UU+ - UU +- 
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if\(u) = 5 2 \ u 



1 + 2u - —II 2 + -u 3 - — u 4 I lnu - Li 2 (u) 
6 4 10 



and 



113 3 

+u [l + 2u u 2 + -u 3 u A I lnu - Li 2 (u) 

6 4 10 



fo ( u ) — $ [ u In u + m 2 In m + uu\ , ip\ (u) = S 2 



23 _ 47 2 _n 41 , , tH 

uu H u m u u H 

60 15 36 6 



m 2 In u + it 2 In U + HU+ — w 2 u 2 



^ 2 ( M )=5 2 



u 2 In u + u 2 In u + uu + —u 2 u 2 



V 2 3 (u) = S 2 



u 2 In u + u 2 In ii 



13 9 9 1 q q 

li U 

12 6 



<p 2 (u)=5 2 



9 , 9 i 2 2 S 3 

u Inu + u inu + uu-\ u u u u 

60 30 



where lA a (x) = I]^ 1 a; n /n a . 

The twist-4 contribution to the FF (|3.7|l can be formulated in terms of the twist-4 specrtal density p^(Q 2 : s): 



F$(Q 2 ) = 



y/2U 



where 



H W iQ 2 )= r^g2 /t ,(4) (Q2s) 
The twist-4 spectral density and the function H^(Q 2 ) have the decompositions 



.,(4) 



(Q 2 , s) = 25 2 {Q 2 ) K n {Q 2 )pi(Q 2 , a) 



n=0 



and 



ff (4) (t? 2 ) = ^E^(Q 2 )^(Q 2 )- 



The explicit expressions for the components of p^(Q 2 , s) and H^(Q 2 ) are written down below (hereafter t = Q 2 ) 

Po(M) 



(s+ty 





s 






is-t + 2sbx 




- 2tln 


_s + t_ 


.*» + *. 



pt(t, S ) = J -— ri \ S 2 -t 2 + 2s 2 ln 



3(s + i) 5 



[6s 3 + 5s 2 t - hst 2 - 6t 3 ] 



s + 1 
2t(s - 2t) 



- 2i 2 In 



In 



s + 1 

2s(2s-t) 
+ (s + ty 



pi(t,s) = 



i 



3(s + i) 

4i(s 3 - ^ 2 - £ 3 ) 
(s + t) 6 



( S + <)4 

y [6s 5 + 19s 4 i + 10s 3 t 2 - 10s 2 t 3 - 19si 4 - 6i 5 ] 
4s(s 3 + s 2 i-i 3 ) 



In 



s + t 



hi 



s + 1 



(s + tf 



■In 



s + 1 
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st(s - t) 
30(7+7) 



Y [21s 2 -65st + 21t 2 



+ — — \t [5s 3 - 5s 2 i - 10si 2 - 6t 3 ] + s [6s 3 + 10s 2 i + 5si 2 - 5t 3 l ) 

3(s + i) b 



2t(3s 3 - s 2 t - st 2 - 2t 3 ) , 
H ; — In 



(s + t) 

Here we introduce new notations 



s + t 



2s(2s 3 + s 2 t + st 2 - 3t 3 ) , 

+ i^TW ln 



s + t 



and get: 



A = \n 



t 



s + 1 



, B = ln 



so 



so +t 



so + i 



— Li 2 , 

6 \s + t 



so 



(so + t) 2 



[t 2 (A - 1) - ts - s (2i + « )B] + S(t), 



Hf(t) = q^^^ [ 12 ^ - 9^ 3 - too(15t + 6s ) - 6s 2 (3t + s )B] + ^ ± ^ + 5(t), 



+ 



6(s + t) 4 9(s + t) 3 (s + i) 3 



[is (s + f)+i 3 A + s 3 B] 



+ T^Ttj 2 ^ ~ 2s °' ) + 2t2A + 2s °^ ~ + 



H 4 3 (t) = 



3(s + i) 6 5(s + t) 5 



{t [3i 4 + 4s (t 3 + 2t 2 s + 2ts 2 + s 3 )] + 4t 5 A + 4s 5 B} 



t A (l + 6A) 3 r4 o „ 22 , , 4 t 3 (-7 + 24A) 

[i 4 + 4i 3 s + 6i 2 s 2 + 4ts 3 + 4s 4 B] 



6(s + i) 4 4(s + t) 4 



9(s + t) 3 



107t 6 t 5 (87+ 100A) 2 r - 4 , 2 , , 4 , -, 

+ — 7^777 ; — + — i — 77? [* + 5i so + I0t 3 sl + lO^s 3 , + 5ts% + 5s 5 Q B] 



90(s + tf 50(s + i) 5 5(s + tf 



t 4 (67+ 150A) 15 r4 



60(s + t) 4 8(s + tf 



[t 4 + 4t 3 s + 6t 2 s 2 + 4ts 3 + 4s^B] + 



3 , 4 „4 D i , i 3 (-49 + 300A) 



90(s + tf 



+ 



40 

9(s +i) 3 



[t 3 + 3i 2 s + 3ts 2 + 3s 3 B] + 2 ^ +t y [*(* - 14s o) + 4t 2 A + 2s (6t - 7s )B] + S{t). 
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